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ABSTRACT 

We report measurements of the cluster X-ray luminosity function out to z = 0.8 based on the final 
sample of 201 galaxy systems from the 160 Square Degree ROSAT Cluster Survey. There is little 
evidence for any measurable change in cluster abundance out to z ~ 0.6 at luminosities less than a 
few times 10 44 ergs s _1 (0.5-2.0 keV). However, between 0.6 < z < 0.8 and at luminosities above 
10 44 

^50 er S s s 1 i the observed volume densities are significantly lower than those of the present-day 
population. We quantify this cluster deficit using integrated number counts and a maximum-likelihood 
analysis of the observed luminosity-redshift distribution fit with a model luminosity function. The 
negative evolution signal is > 3cr regardless of the adopted local luminosity function or cosmological 
framework. Our results and those from several other surveys independently confirm the presence 
of evolution. Whereas the bulk of the cluster population does not evolve, the most luminous and 
presumably most massive structures evolve appreciably between z = 0.8 and the present. Interpreted 
in the context of hierarchical structure formation, we are probing sufficiently large mass aggregations 
at sufficiently early times in cosmological history where the Universe has yet to assemble these clusters 
to present-day volume densities. 

Subject headings: cosmology: observations — galaxies: clusters: general — X-rays: general 



1. INTRODUCTION 

Galaxy clusters form via the gravitational amplifica- 
tion of rare, high peaks in the cosmic matter density field. 
The redshift evolution of cluster abundance depends on 
the growth rate of density perturbations which is, in turn, 
sensitive to the mean cosmic matter density (Ojw) and, 
to a lesser extent, the dark-energy density (Oa)- Thus 
observations of cluster space density with sufficient tern- 



straints (e.g., 


Oukbir & Blanchard 1992: Eke et alll998t 


Bahcall ct al. 


1999). 



A directly measurable and robust diagnostic of cluster 
abundance is the X-ray luminosity function (XLF) , that 
is the volume density of clusters per luminosity interval. 
X-ray selected galaxy clusters are particularly well-suited 
for this type of analysis. Clusters are efficiently detected 
at X-ray wavelengths to high redshift (currently out to 
z ~ 1.25) thus provi ding the leverage for evolutionar y 
studies (see review bv lRosati. Borgani. fc Normanl 2002). 
The resulting samples feature high statistical complete- 
ness which is clearly important for deriving reliable num- 
ber counts. Since X-ray surveys have well-defined se- 
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lection functions, it is straight-forward to convert these 
number counts into volume-normalized measures such as 
the luminosity function. Finally, given the strong corre- 
lation between X-ray luminosity and cluster mass, it is 
possible to transform the observed XLF into the cluster 
mass function which is the fundamental relation in the 
theoretical treatment. 

Though the cluster XLF was first measured with 
an X-ray flux-limited sample over two decades ago 
l)Piccinotti et al J 1198 2). a definitive characterization of 
its evolution has proven very difficult. The latter is a 
particularly important test in observational cosmology 
because the variation of the cluster XLF as a function 
of redshift reflects the evolution of the cluster mass func- 
tion. Such a measure allows for strong constraints on Qm 
even account ing for uncertainty in the luminosity-mass 
relation fe-g.-IBorgan i et all 2001). Early theoretical pre- 
dictions (e.g., Kaiser 1986) postulated dramatic positive 
evolution in the XLF where the volume density of clus- 
ters of a fixed luminosity would increase with redshift. 
This would be an "observer-friendly" universe since the 
loss in sensitivity at high redshift in a flux-limited survey 
would be offset by the growing population of detectable 
sources. Contrary to this scenario, observations of the 
cluster XLF range from zero evolution to negative evolu- 
tion depending on the redshifts and luminosities probed. 
These findings are consistent with current predictions for 
a low-density universe where mild negative evolution is 
restricted to the most luminous, high-redshift clusters 
while there is little change in the bulk of the population 
fe.g.. lBorgani fc Guzzdl200ll and references therein). 

Taking advantage of the first time an X-ray sur- 
vey extended into cos mologically interest ing redshifts, 
iCioia et all l|1990al) and lHenrv et all (|1992D used 67 clus- 
ters (0.14 < z < 0.60) from the pioneering Einstein Ex- 
tended Medium Sensitivity Survey fEMSS: ICioia et alJ 
MM IStncke et al.l [EM IMaccacaro et alJ 119941) to 
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Fig. 1. — X-ray lumin osity and red shift distributi on of the 160SD ijMullis et al.H2003) . the EMSS IGioia fc LuDpinolfl994l with updates 
from the literature, e.g.. IHenrv 2000), and the BCS lEbclhiH ct al. 1998) cluster samples. The dotted curves (left to right) are indicative 
flux limits of 2.7 X 10~ 12 , 1.5 X 10~ 13 , and 3 X 10~ 14 ergs cm " 2 s" 1 (0.5-2.0 keV) The EMSS and BCS luminosities were converted to 
this energy band assuming a Raymond-Smith plasma spectrum (Raymond & Smith 1977) with a metallicity of 0.3 solar and the reported 
gas temperature (either directly measured or estimated from the luminosity-temperature relation). 



make the first detection of evolution in the observed 
XLF. Based on a steepening of the high-redshift lumi- 
nosity function, they found a deficit of high-luminosity 
clusters at z > 0.3 with a statistical significance of ap- 
proximately 3d. Though most subsequent investigations 
corroborate these findings, some questions have been 
raised concerni ng the reliability of the EMSS evolution 
detection (e.g.. INichol et alJll997t lEbehng et alJl2000bl 
Ellis & Jonesl2002HLewis et al.l2002j) . Of historical note, 
the only other pre- ROSA T measurement of significant 
cluster evolution came from lEdge et alJ l)1990f) who found 
rapid negative evolution in the luminosity function at 
z < 0.2 based on a HEAO-1 sample. This was ultimately 
overruled b y a definitive and no n-evolving measure of the 
local XLF (jEbehng et aTIIT997|) . 

Seeking in part to confirm or refute the EMSS's contro- 
versial claim of negative cluster evolution, a large num- 
ber of cluster sur veys were in i tiated in the 1 990s based 
on ROSAT data <IVogeslll992t ITmmperllT993l) . The 160 
Square Degree ROSAT Cluster Survey (h ereafter 160SD, 
Vikh linin et alJl998atlMullis et al.12 003) is one such pro- 
gram and the subject of this paper. Additional sur- 



veys probing to high re dshift include the NEP (Mullis 
l200ll IHenrv et al.N2001t IGioia et al l 12001 BMW-HRI 



lIMoretti et all 120011: IPanzera et al J 1200*31. BSHARC 
iRomer et alJ 120001). MACS iEbeling Edge fc Henrvl 
1200111 . RDCS URosati et"aDll99flt 119981 120001) . RIXOS 
llCastander et alJ 1199.1 IMason et al.l "*2000). SSHARC 



iBurke et all l200l. and WARPS l|Scharf et all Il997t 
IPerlman et al.ll2002j) . Complementary surveys at low 
redshifts (z < 0.3) yielded accurate determinations of 
the local luminosity function, thus providing the crucial 
baseline for the detection of redshift evolution. These lo- 
cal sur veys include t he BCS+eBCS llEbeling , et alJH99Sl 
2000a), RASS1BS iDe Grandi et all Il999al). and RE- 



FLEX ( Bohrin ger et alJ200Hl . From the shear number of 
projects it should be clear that ROSAT was a watershed 
event for X-ray cluster surveys. 

Critical comparisons of the ROSAT and EMSS lumi- 
nosity functions must account for the overlap (or lack 
thereof) of the measurements in redshift and luminos- 
ity. Insufficient attention to this point led to confusion 
in some early analyses and debates. Different combina- 
tions of sensitivity and areal coverage in flux-limited sur- 
veys, when convolved with the intrinsic cluster luminos- 
ity function, result in populating different regions of the 
observed luminosity-redshift plane (see indicative results 
in Figure QJ. In practice it is difficult to directly test the 
cluster evolution seen in the EMSS detection because it 
lies at the extreme bright end of the luminosity function. 
Thus large search volumes are required to detect ade- 
quate numbers of such rare clusters. The 160SD is one 
of the largest serendipitous X-ray surveys conducted with 
ROSAT. With this substantial areal coverage and high 
sensitivity, the 160SD survey is well-positioned to probe 
cluster evolution. Preliminary analyses of our sample, 
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flux (0.5-2.0 keV) [ergs s~ 1 cm' 2 ] 



Fig. 2. — The selection function of the 160SD survey specifying 
the square degrees of area covered as a function o f limiting X-ray 
flux (0 .5-2.0 keV). These data are tabulated in IVikhlinin et ahl 
(1998a, Table 5). The shading encodes this information for subse- 
quent use in Figure |3] The vertical line indicates the minimum flux 
limit (4 X 10 4 ergs cm -2 s —1 , 0.5— 2.0keV) used in the present 
analysis. 



with at-the-time incomplete optical follow-up, seemed 
to confirm the deficit of high-lu minosity, high-r edshift 
cluste rs first seen by the EMSS ijVikhlinin et alJll998rl 
2000). 

In this paper we present measurements of the clus- 
ter XLF out to z = 0.8 and describe their implications 
for cluster evolution based on the final 160SD sample of 
201 clusters. In § 2 we outline the basic construction of 
our cluster sample and the associated selection function 
used herein. The formalism associated with the XLF is 
defined in § 3 and used to measure cluster abundances 
at both low and high redshifts. In § 4 we characterize 
the evolution in the cluster population using integrated 
number counts and a maximum- likelihood analysis of the 
observed luminosity-rcdshift distribution relative to an 
evolving Schechter function. We discuss our results in 
the context of previous works and draw conclusions in § 5. 
Throughout this analysis we use the cosmological param- 
eters H = 50 h 50 km s _1 Mpc -1 , Qm — L an d = 
(Einstein-de-Sitter, EdS) for direct comparison to pre- 
vious work in this field. We also repeat calculations in 
the currently preferred cosmology where J7m = 0.3 and 
Q\ = 0.7. X-ray fluxes and restframe luminosities are 
quoted in the 0.5-2.0 keV energy band unless otherwise 
stated. 

2. THE 160SD CLUSTER SAMPLE 

The 160SD sample of 201 galaxy clusters is the largest 
high-redshift, X-ray selected sample published so far. For 
instance, there are 73 objects at z > 0.3 and 22 objects 
at z > 0.5. The 160SD sample was constructed via the 
serendipitous detection of extended X-ray sources in 647 
archival ROSAT PSPC observations. Of 223 cluster can- 
didates, we identified 201 as galaxy clusters, 21 as prob- 
able false-detections due to blends of unresolved point 
sources, and one source is unidentified due to its proxim- 
ity to a bright star. We have se cured spectroscopic red- 
shifts for 200 of the 201 clusters. IVikhlinin et alJ (|1998a^ 
give a complete descri ption of the survey methodol- 
ogy. iMullis et al.l l)2003|) detail the optical follow-up and 
present the final cluster catalog with spectroscopic red- 



shifts. Note that the number of false-detections in the 
sample agrees very well with that expected from the con- 
fusion of point sour ces as demonstrated in the Monte- 
Carlo simulations of IVikhlinin et all l)1998a|L Two false- 
detections (Nos. 77 and 141) serendipitously imaged by 
XMM-Newton and Chandra, respectively, are indeed con- 
fused point sources. Finally, the optical survey imaging 
was sufficiently deep to demonstrate that none of the 
false-detections should be galaxy clusters at z < 0.9. 

Although spatial X-ray extent was our primary selec- 
tion criterion, detailed comparisons with other surveys 
demonstrate that no known clusters were miss ed as un- 
resolved sources (see §4 of lMnllis et aTl 120031 and ref- 
erences therein). Thus the 160SD clusters are in ef- 
fect drawn from a statistically complete, flux-limited 
survey with an areal coverage (fi) and sensitivity char- 
acterized by the selection function shown in Figure 
A total of 158.5 deg 2 were surveyed and the me- 
dian survey flux (where = Sltotai/2) is 1.2 x 10~ 13 
ergs cm -2 s -1 (0.5-2.0keV). We restrict our analysis to 
a minimum flux of 4 x 10~ 14 ergs cm" 2 s" 1 to avoid 
any uncertainties in the sky coverage at very faint fluxes. 
There are 190 clusters in the 160SD survey sample above 
this flux limit. Figure ^shows the position in luminosity- 
rcdshift space of our sample relative to the BCS, one of 
the key reference samples at low redshift, and the EMSS. 
Figure shows the 160SD data in greater detail. Here 
the cluster luminosities are plotted as a linear function 
of comoving volume which provides a more uniform view 
of the volume sampling. The greyscale in the figure indi- 
cates the parameter space probed by our survey — sen- 
sitivity is maximal in light regions and minimal in dark 
regions. 

In our subsequent derivations of the cluster XLF and 
tests for evolution, we exclude several objects (all at z < 
0.3) to minimize the biasing of our results. To avoid 
potentially skewing the impartiality of the sampling, we 
reject nine clusters whose redshifts are within Az = 0.015 
of the original target of the ROSAT PSPC observations 
(Nos. 16, 32, 112, 134, 165, 166, 174, 177, 206). Four 
X-ray-overluminous elliptical galaxies or "fossil groups" 
(Nos. 110, 144, 201, and 211) are discounted because 
these special structures are unlikely to meet the se lection 
criteria of the local samples (|Vikhlinin et alll999lL Thus 
the final cluster sample used here consists of 177 clusters. 

3. THE X-RAY LUMINOSITY FUNCTION 

We define the cluster differential luminosity function 
to be 



<p(L x ,z) 



d 2 N 
dVdLx 



(1) 



where N is the number of clusters of luminosity Lx in 
a volume V at a redshift z. The standard approach 
for deriving a nonparametric representation of the dif- 
ferential cluster X LF is bas e d on the 1/V max technique 
first proposed by iSchmidtl l)1968l) and generalized by 
lAvni fc Bahcalll l)1980jh Here the observed luminosity 
range is parsed into bins of width AL, within each of 
which there are Nj observed clusters. The XLF is es- 
timated by summing the density contributions of each 
cluster in the considered luminosity/redshift bin, 
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Fig. 3. — X-ray luminosity of the 160SD clusters versus redshift z (top axis) and comoving volume out to redshift z (bottom axis). 
The volume corresponds to a survey solid angle of 158.5 deg 2 in an Einstein— de-Sitter cosmology. The shading, established in Figure l2l 
demonstrates the range of survey flux limits and associated sky coverage in this parameter space. The thick, solid curve indicates the 
minimum flux limit (4 X 10~ 14 ergs cm~ 2 s -1 , 0.5-2.0keV) used in the present analysis. The horizontal (vertical) dashed lines are the 
boundaries of the luminosity (redshift) bins used in computing the luminosity function. 



1 j 1 

where V max is the total comoving volume in which a clus- 
ter of luminosity Lx,i could have been detected above the 
flux limits of the survey. Over a specific redshift interval 
(^min < z < z max ), this search volume is defined to be 

= r* x n(f x (L XlZ ))^±dz. (3) 

Here 0(/x) is the sky area surveyed in steradians as a 

function of X-ray flux and dV } is the differential, co- 
moving volume element p er steradian. 

iPage fc Carreral l|200fD describe a refinement to the 
canonical approach in which the estimator takes the form 

4>(Lx,z) = f L X max fZmax . \ \ dV(z) j , r 

. W 

where the boundaries of the luminosity bin are specified 
by Lx.min and ix.max- The effect of bringing the lumi- 
nosity interval into the double integral results in a better 
estimate of the effective AL which can be smaller than 
the full bin width for regions of luminosity-redshift space 



transected by the survey flux limit (e.g., faint luminosi- 
ties). We will use the Page-Carrera (PC) estimator in our 
subsequent derivations of the cluster XLF. However, the 
results are essentially identical to the classical 1/Vmax 
procedure except at the very faint end where the XLF is 
marginally increased. 

As for a parametric representation of the cluster 
XLF, observatio ns are well fit by a Schechter function 
ljSchechterlll97fl) of the form 

HLx,z)dL x = ^)" Q exp(-^)^ (5) 

where <j>* is the normalization (units h\ Mpc~ 3 ), a is 
the faint-end slope, and L x is the characteristic luminos- 
ity marking the interface between the power-law and the 
exponential regimes. An equivalent expression for the 
XLF commonly used in the literature is 

4>{Lx,z) = AL x a exp(^-^j . (6) 

With Lx in units of 10 44 ergs s _1 , the as- 

sociated normalization, A, has units of /ifg Mpc -3 
(10 44 h^Q ergs s _1 ) Q_1 , and the two normalizations are 
related by <f>* = A{L\f- a . 
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Table 1. Best-Fit Schechter Parameters for the Local Cluster XLF 
(z < 0.3) 



Sample 




a 


<t>* 


Reference 




(10 44 ergs s" 1 ) 




h\ Q Mpc -3 





BCS 5.70±J || 1.85±{5;g! (7.56± ;?5) x 10" 8 Ebeling et al. (19971 

RASS1BS 3.80±g;^ 1.52±g;^ (2.53± ;^) x 10" 7 De Grand! et al. (1999b~) 

REFLEX* 4.21+n?4 1.63j;{5;2g C 1 - 80 ^'.! ) x 10 ~ 7 Bohringer et al. (20021 

5 - 18 ±o'.50 L69 iao45 lm x 10" 7 Q M = 0.3, fi A = 0.7 



Note. — Parameters for an Einstein— de-Sitter universe unless otherwise indicated 
t_L^ is quoted in the 0.5-2.0 keV band 

^ fit includes a correction for missing flux for consistency with the other surveys (see details in Bohringer et al. 2002) 



3.1. The Local XLF 

Knowledge of the local or near present-day abun- 
dance of clusters is fundamental to evolutionary stud- 
ies because it serves as the no-evolution baseline against 
which distant cluster samples can be tested. As previ- 
ously noted, one of the significant achievements of the 
ROSAT era is the accurate determination of the lo- 
cal cluster XLF (z < 0.3). Three principal measure- 
ments are based on the BCS, RASS1BS, and REFLEX 
samples which were constructed by surveying large por- 
tions of the two-thirds of the sky outside the Galactic 
plane 9 at relatively bright fluxes. lEbeling et al.1 l|1997fl 
reported the first results based on 199 BCS clusters in the 
northern hemisphere (fx > 2.8 x 10 -12 ergs cm -2 s -1 , 
0.1-2.4 keV). P art of a pilot pr o gram o f the larger 
REFLEX survey, iDe Grandi et alJ l)1999bj) presented a 
measurement based on 129 RASS1BS clusters (fx > 
3-4 x 10 -12 ergs cm -2 s -1 , 0.5-2.0keV) from the south 
Galactic cap. Finally, the XLF determination based on 
the l argest sample to date c omes from the REFLEX sur- 
vey. iBohringcr et alJ l|2002fl made a detailed analysis of 
452 clusters extracted from the southern celestial hemi- 
sphere (fx > 3 x 10~ 12 ergs cm~ 2 s _1 , 0.1-2.4keV). 

Nonparametric determinations of the local XLF (z < 
0.3) from the all-sky samples (BCS, RASS1BS, and RE- 
FLEX) are shown in Figure We also plot the best- 
fitting Schechter functions for these data and list the 
associated best-fit parameters in Tabled for future ref- 
erence. The results demonstrate that we have accurate 
knowledge of the local cluster luminosity function. Inde- 
pendent investigators using different X-ray selection pro- 
cedures over different regions of the sky agree on the local 
space density of clusters. For example, in the luminos- 
ity interval 5 x 10 43 - 10 45 ergs s" 1 (0.5-2.0 keV), the 
internal accuracy of these XLF measurements is approx- 
imately ±10%-20% (estimated from the ±lcr excursion 
of the error envelopes plotted in Figure0J . Moreover, the 
systematics are also small — the results from the BCS 
and RASS1BS surveys vary a maximum of about ±25% 
relative to the Schechter fit of REFLEX. 

9 Note that the significant gap in all-sky coverage due to the 
former zone of avoidance (|fe| < 20°) is be ing redressed by the 
CIZA survey lEbeling. Mullis. fc TulTvll2'002Tl 



Before considering the cluster population at high red- 
shift, we first examine the low-redshift diagnostic using 
the 160SD cluster sample. Using the procedure pre- 
viously described, we estimate the local XLF between 
~10 42 and 3 x 10 44 ergs s~ 4 using the 110 clusters at 
0.02 < z < 0.3 in the 160SD survey. Our measurement 
is plotted in Figure^Jand is tabulated in Tabled The lu- 
minosity binning is uniform in log space and data points 
are plotted at the center of the luminosity interval. The 
error bars are the equivalent ±ler uncertainties based 
on Poissonian errors for the number of clusters per bin 
l)OehrelsE p86). Note the good agreement at low redshift 
between the 160SD and the three principal local samples. 
In Figure [5] we add the measurements from four addi- 
tional deep surveys (RDCS, EMSS, NEP, and WARPS) 
thus producing a compilation of all of the ROSAT (plus 
EMSS) local XLFs published to date from X-ray selected, 
X-ray flux-limited cluster samples. The shallow all-sky 
surveys accurately measure the local luminosity function 
at intermediate to high luminosities (> 10 43 ergs s _1 ) 
but are relatively insensitive to very low luminosity clus- 
ters (< 10 43 ergs s _1 ). In a complementary fashion, the 
deep surveys better measure the faint end of the local 
XLF, provide reasonable precision at intermediate lumi- 
nosities, but poorly constrain the very bright end (> 10 44 
ergs s _1 ) due to limited survey volumes at low redshift. 
The cluster luminosity-redshift distributions in Figure ^ 
also illustrate this dependence on flux limit and solid an- 
gle. 

3.2. The High-Redshift XLF 

We measure the distant XLF using the 66 clusters 
from the 160SD sample at 0.3 < z < 0.8 and fluxes 
above 4 x 10~ 14 ergs cm~ 2 s" 1 . The numerical results 
are shown in Table Eland plotted in Figure El We have 
probed a sufficiently large volume such that we can de- 
rive useful results in two intervals: 0.3 < z < 0.6 and 
0.6 < z < 0.8. In addition to our high-redshift measure- 
ments, we also show in FigureHOthc local determinations 
of the XLF. These establish the regime where the high- 
redshift results should lie if the spatial density of clusters 
does not evolve out to the considered redshifts. 

Our measurement of the XLF at 0.3 < z < 0.6 
probes the luminosity range 2 x 10 43 - 7 x 10 44 
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L x (0.5-2.0 keV) [ergs s" 1 ] 

Fig. 4. — Determinations of the local cluster X-ray luminosity function as measured by the 160SD survey and the local reference samples 
(BCS, RASS1BS, and REFLEX) in an Einstein-de-Sitter universe. The 160SD data values along with the number of clusters and average 
clust er redshift for each luminosity bin are given in Table El Nonpa rametric d ata points and Schechter fits for the reference samples are 
from lEbeling et al.l tl99^ . IDe Grandi et alJ I1999bl) . and lBr3nringer et alj|2002TV respectively. The BCS data points arc based on a merged 
analysis of the BCS+eBCS samples (H. Ebeling 2003, private communication). The shaded region indicates the ler uncertainty envelope 
of the Schechter fits assuming the errors on L^r and a are correlated. The indicated uncertainties on data points are ±1ct. 




Fig. 5. — Compilation of local XLFs as mea sured by eight X-ray flux-limited surveys. RDCS: Rosati et al. (1998), EMSS: Henr v^t alJ 
tl992T) . NEP: lGloia et alJ Pop!!) , and WARPS: Pones et alJ 12000ri) and the references in Figure 4 (Einstein-de-Sitter universe). 



CLUSTER XLF EVOLUTION 



7 



T 1 1 1 1 1 1 1 — | 1 1 1 1 1 1 1 1 — |" 




L x (0.5-2.0 keV) [h^ 2 ergs s" 1 ] 

Fig. 6. — The high-redshift cluster X-ray luminosity function from the 160SD sample (Einstein-de-Sitter universe). The data points for 
0.6 < z < 0.8 have been slightly offset to avoid confusion. The 160SD data values along with the number of clusters and average cluster 
redshift for each luminosity bin are given in Table [3] Schechter fits to the local XLF are also plotted as described in Figure[I] If no clusters 
are detected in a luminosity bin, we plot the 1<t upper limit for the Poisson error on zero clusters. 




1Q 43 1Q 44 1Q 45 

L x (0.5-2.0 keV) [h^ 2 ergs s" 1 ] 

Fig. 7. — The high-redshift cluster X-ray luminosity function from the 160SD sample for a cosmology characterized by the parameters 
Qm = 0.3 and = 0.7. The data points for 0.6 < z < 0.8 have been slightly offset to avoid confusion. The 160SD data values along 
with the number of clusters and average cluster redshif t fo r each luminosity bin are given in Table El The REFLEX XLF was explicitly 
measured for a A-dominated cosmology by Bohringcr ct al. (2002). We transform the BCS and RASS1BS results to this cosmology based 
on the REFLEX results (see text for details). 
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Table 2. Local 160SD Cluster X-ray Luminosity Function (0.02 < z < 
0.3, EdS) 



Lx(min) t L x (max) t <t>(Lx) 4>( L x)-l<r 4>( L x)+la < z > N c] 

(fc| Mpc~ 3 (10 44 ergs s^ 1 )- 1 ) 



Lx (center) ' 
(10 44 ergs s' 1 ) 



0.016 


0.013 


0.019 


5.09 


X 


10~ 4 


1.78 


X 


10~ 4 


1.18 


X 


10~ 3 


0.045 


2 


0.022 


0.019 


0.026 


2.19 


X 


10~ 4 


7.65 


X 


10~ 5 


5.07 


X 


10~ 4 


0.035 


2 


0.031 


0.026 


0.037 


1.42 


X 


io~ 4 


6.45 


X 


10~ 5 


2.81 


X 


10~ 4 


0.062 


3 


0.044 


0.037 


0.053 


2.08 


X 


10 


3.59 


X 


10 -6 


6.85 


X 


10 ° 


0.135 


1 


0.063 


0.053 


0.075 


4.59 


X 


io- 5 


2.61 


X 


IO" 5 


7.69 


X 


IO" 5 


0.142 


5 


0.089 


0.075 


0.106 


2.87 


X 


io- 5 


1.80 


X 


IO" 5 


4.42 


X 


IO" 5 


0.136 


7 


0.125 


0.106 


0.149 


1.86 


X 


10~ 5 


1.28 


X 


10~ 5 


2.65 


X 


IO" 5 


0.156 


10 


0.177 


0.149 


0.211 


6.98 


X 


10~ 6 


4.55 


X 


10~ 6 


1.04 


X 


10~ 5 


0.154 


8 


0.251 


0.211 


0.298 


7.44 


X 


10~ 6 


5.65 


X 


10~ 6 


9.72 


X 


10~ 6 


0.200 


17 


0.355 


0.298 


0.422 


5.49 


X 


10~ 6 


4.35 


X 


IO" 6 


6.89 


X 


IO" 6 


0.197 


23 


0.502 


0.422 


0.597 


1.27 


X 


10~ 6 


8.55 


X 


IO" 7 


1.86 


X 


IO" 6 


0.195 


9 


0.710 


0.597 


0.844 


9.89 


X 


io- 7 


6.94 


X 


io- 7 


1.39 


X 


IO" 6 


0.239 


11 


1.004 


0.844 


1.194 


3.60 


X 


io- 7 


2.16 


X 


io- 7 


5.76 


X 


io- 7 


0.226 


6 


1.420 


1.194 


1.688 


8.26 


X 


io- 8 


2.89 


X 


ur 8 


1.92 


X 


IO" 7 


0.225 


2 


2.008 


1.688 


2.388 


5.76 


X 


IO" 8 


2.01 


X 


IO" 8 


1.34 


X 


IO" 7 


0.180 


2 


2.840 


2.388 


3.377 


2.01 


X 


IO" 8 


3.47 


X 


IO" 9 


6.63 


X 


IO" 8 


0.296 


1 



tix i s quoted in the 0.5-2.0 keV band 



ergs s _1 (Figure El filled circles). Except for the data 
point near 5 x IO 43 ergs s _1 which is 1.9a off the local 
relation, there is excellent agreement between these 
results and the no-evolution benchmark at least out to 
intermediate luminosities, Lx ~2 x IO 44 ergs s _1 . In 
this region the 160SD best matches the normalization of 
the REFLEX XLF. At higher luminosities the distant 
XLF is lower than the local population but we will 
demonstrate in S I4.1I that this effect is only marginally 
significant with respect to the local XLF with the 
highest normalization (RASS1BS). The median redshift 
for these depressed data points is < z >= 0.50 (see 
Table El . 

At the highest rcdshifts probed by the 160SD in 
the present analysis, 0.6 < z < 0.8, we measure use- 
ful constraints over the luminosity interval 1 - 6 x 10 44 
ergs s _1 , and the median cluster redshift is < z >= 0.70. 
(Figure El open squares). The distant cluster volume 
densities are below the local level at all measured lu- 
minosities. We will show that this result is significant 
(> 3cr) in both the Einstein-de-Sitter and A-dominated 
models even with respect to the lowest-normalization lo- 
cal XLF (§ m . This is a deficit of high luminosity, 
high redshift clusters in a manner similar to that seen 
in the EMSS results. Note that although the optical 
CCD imaging used to identify the 160SD clusters was 
sufficiently deep to reveal massive clusters to z = 0.9, we 
have conservatively cut off the most distant redshift shell 
at z = 0.8. This lessens any potential negative evolution 
in the results since the larger redshift boundary would 
increase the search volume without adding clusters and 
hence depress the data points. 

To assess the impact of changing the cosmological 
framework from an Einstein-de-Sitter to a A-dominated 
universe, we repeat our calculation of the cluster lu- 
minosity function setting Qm = 0.3 and = 0.7, 
the results of which are shown in Figure [7| and Table 0] 
The Einstein-de-Sitter XLF (4>EdS, Figure EJ) and A- 
dominated XLF (</>a) appear very similar because the 



data points of the latter are offset diagonally down and 
to the right approximately along the former. This is 
the combined effect of the increase in both the cluster 
luminosities and the search volumes in a A-dominated 
universe. The actual positioning of the two XLFs rela- 
tive to each other depends on the cluster redshifts and 
the specific luminosity interval. For example in the RE- 
FLEX survey (< z >= 0.08), between 10 43 and 10 45 
ergs s _1 the ratio of the fitted XLFs (</>A/</>Eds) is less 
than unity with a broad minimum of ^0.8 around 10 44 
ergs s~ x . REFLEX is the only local sample for which 
the XLF in a A-cosmo logy has been explicitly measured 
(Boh ringer et alJ l2002) . However, given the similarity of 
the redshift and luminosity distributions, we have used 
the ratio </>A/<^EdS from REFLEX to make an approxi- 
mate transform of the BCS and RASS1BS to this alter- 
nate cosmology. 

It is apparent from Figured that the high-redshift 
160SD XLFs and the local XLFs shift in similar ways; 
thus the apparent deficit of high-luminosity clusters per- 
sists in the A-dominated cosmology. The only im- 
portant difference, of course, is that the point where 
the 160SD data depart significantly from the non- 
evolution baseline is shifted to larger luminosities, and 
this is anticipated given the increase in luminosity dis- 
tance (L x > 3 x 10 44 ergs s" 1 for 0.3 < z < 0.6, and 
L x > 2 x 10 44 ergs s" 1 for 0.6 < z < 0.8). 

We will examine the significance and strength of this 
apparent cluster evolution in the following section. 

4. QUANTIFYING EVOLUTION 

X-ray luminosity functions like those shown in 
Figures 0H3 are useful for visualizations and qualitative 
assessments of the cluster population; however, they are 
non-optimal for quantitative analyses. For example am- 
biguities exist in the selection of the luminosity binning 
(e.g., fixed or adaptive intervals) and the loci of the plot- 
ted data points in luminosity space (e.g., at the bin center 
or the density- weighted mean luminosity). Furthermore, 
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Table 3. Distant 160SD Cluster X-ray Luminosity Function (EdS) 
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Lx (center) t 
(10 44 ergs s- 



Lx(min) t L x (max) t 



(ft| Mpc 



' 3 (10 44 ergs s" 



<f>(Lx)+i<r <z> N ci 



0.3 < z < 0.6 




























0.251 


0.211 


0.298 


2.86 


X 


10~ 6 


4.95 


X 


io- 7 


9.45 


X 


10~ 6 


0.329 


1 


0.355 


0.298 


0.422 


2.01 


X 


10~ 6 


9.10 


X 


io- 7 


3.96 


X 


10~ 6 


0.334 


3 


0.502 


0.422 


0.597 


2.80 


X 


10~ 6 


2.06 


X 


10~ 6 


3.77 


X 


10~ 6 


0.376 


14 


0.710 


0.597 


0.844 


5.83 


X 


io- 7 


3.80 


X 


io- 7 


8.72 


X 


io- 7 


0.420 


8 


1.004 


0.844 


1.194 


3.17 


X 


io- 7 


2.18 


X 


io- 7 


4.53 


X 


io- 7 


0.461 


10 


1.420 


1.194 


1.688 


1.61 


X 


io- 7 


1.10 


X 


io- 7 


2.29 


X 


io- 7 


0.418 


10 


2.008 


1.688 


2.388 


9.18 


X 


io- 8 


6.32 


X 


io- 8 


1.31 


X 


io- 7 


0.517 


10 


2.840 


2.388 


3.377 


2.88 


X 


HT 8 


1.64 


X 


IO" 8 


4.83 


X 


io- 8 


0.501 


5 


4.016 


3.377 


4.776 


7.67 


X 


HT 9 


2.68 


X 


io- 9 


1.78 


X 


io- 8 


0.485 


2 


5.679 


4.776 


6.754 


4.85 


X 


10- 9 








4.85 


X 


io- 9 







0.6 < z < 0.8 




























1.004 


0.844 


1.194 


3.65 


X 


io- 7 








3.65 


X 


io- 7 







1.420 


1.194 


1.688 


5.17 


X 


IO" 8 


8.95 


X 


io- 9 


1.71 


X 


io- 7 


0.699 


1 


2.008 


1.688 


2.388 


1.88 


X 


io- 8 


3.25 


X 


io- 9 


6.20 


X 


nr 8 


0.625 


1 


2.840 


2.388 


3.377 


1.59 


X 


IO" 8 








1.59 


X 


io- 8 







4.016 


3.377 


4.776 


4.71 


X 


io- 9 


8.15 


X 


io- 10 


1.55 


X 


io- 8 


0.700 


1 


5.679 


4.776 


6.754 


5.33 


X 


10- 9 








5.33 


X 


10- 9 








f Lx is quoted in the 0.5-2.0 keV band 



Table 4. Distant 160SD Cluster X-ray Luminosity Function (SIm = 0.3 
and Qa = 0.7) 



Lx(min) t L x (max) t <t>(L x ) </>(L x )-i CT <j>(L x ) +la < z > N cl 

(h| Mpc" 3 (IO 44 ergs s- 1 )- 1 ) 



Lx (center) t 
(IO 44 h M 2 ergs s- 1 ) 



0.3 < z < 0.6 
























0.355 


0.298 


0.422 


9.52 


X 


io- 7 


1.65 x 10~ 7 


3.14 


X 


10~ 6 


0.329 


1 


0.502 


0.422 


0.597 


9.51 


X 


io- 7 


4.97 x 10~ 7 


1.70 


X 


10~ 6 


0.351 


4 


0.710 


0.597 


0.844 


1.04 


X 


10~ 6 


7.62 x 10~ 7 


1.40 


X 


10~ 6 


0.372 


14 


1.004 


0.844 


1.194 


1.93 


X 


io- 7 


1.21 x 10~ 7 


2.97 


X 


io- 7 


0.429 


7 


1.420 


1.194 


1.688 


1.20 


X 


io- 7 


8.22 x 10~ 8 


1.71 


X 


io- 7 


0.439 


10 


2.008 


1.688 


2.388 


5.99 


X 


io- 8 


4.12 x 10~ 8 


8.56 


X 


io- 8 


0.440 


10 


2.840 


2.388 


3.377 


3.05 


X 


io- 8 


2.05 x 10~ 8 


4.45 


X 


io- 8 


0.509 


9 


4.016 


3.377 


4.776 


1.47 


X 


io- 8 


9.26 x 10~ 9 


2.27 


X 


io- 8 


0.507 


7 


5.679 


4.776 


6.754 


1.39 


X 


io- 9 


2.41 x 10- 1 " 


4.59 


X 


io- 9 


0.516 


1 


8.032 


6.754 


9.552 


1.75 


X 


10- 9 




1.75 


X 


io- 9 







0.6 < z < 0.8 
























1.420 


1.194 


1.688 


1.59 


X 


io- 7 




1.59 


X 


io- 7 







2.008 


1.688 


2.388 


3.75 


X 


io- 8 




3.75 


X 


io- 8 







2.840 


2.388 


3.377 


6.76 


X 


10- 9 


1.17 x 10~ 9 


2.23 


X 


io- 8 


0.699 


1 


4.016 


3.377 


4.776 


2.90 


X 


io- 9 


5.02 x IO" 10 


9.57 


X 


io- 9 


0.625 


1 


5.679 


4.776 


6.754 


1.51 


X 


io- 9 


2.61 x IO" 10 


4.98 


X 


io- 9 


0.700 


1 


8.032 


6.754 


9.552 


1.65 


X 


10- 9 




1.65 


X 


10- 9 








tLx is quoted in the 0.5-2.0 keV band 



for the case of negative evolution the effect that we are 
attempting to measure is either a diminishing signal or 
a non-detection. Thus we will apply the alternate ap- 
proaches of integrated number counts and a maximum- 
likelihood fit of an evolving model XLF to quantify the 
significance and strength of the apparent evolution in the 
160SD clusters. 



4.1. Integrated Number Counts 

For a given region of luminosity-redshift space we com- 
pare the number of observed clusters {N ohs ) with the 
number that are expected (iV eX p) assuming there is no 
evolution in the population. The latter is computed by 
integrating the local luminosity function, <fi(Lx,z), over 
luminosity and redshift, and folding this through the sur- 
vey selection function, 0(/x), using the following equa- 
tion, 
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Table 5. Observed versus Expected Number of Clusters (L x , m i n < 
L x < oo) 





N ohs 


Wcxp 


significance (ct) 


(10 44 h-£ ergs s" 1 ) 




REFLEX BCS RASS1BS 


REFLEX BCS RASS1BS 



0.3 < z < 0.6, EdS 



4.776 





1.5 


3.4 


5.2 


2.3 


1.8 


2.6 


3.377 


2 


9.1 


6.7 


10.8 


2.5 


1.8 


3.0 


2.388 


7 


16.2 


11.8 


19.2 


2.1 


1.3 


3.0 


1.688 


17 


25.8 


18.7 


30.4 


1.7 


0.2 


2.5 


< 0.6, n M = 


0.3 and 


= 0.7 












6.754 





3.5 


3.0 


4.0 


1.9 


1.6 


2.1 


4.776 


1 


7.8 


6.1 


9.1 


2.7 


2.2 


3.1 


3.377 


8 


14.8 


11.0 


17.4 


1.7 


0.7 


2.3 


2.388 


17 


24.4 


17.8 


28.8 


1.4 


0.0 


2.2 


< 0.8, EdS 
















4.776 





4.1 


3.2 


4.8 


2.1 


1.7 


2,1 


3.377 


1 


7.9 


5.8 


9.3 


2.7 


2.0 


3.1 


2.388 


1 


12.6 


9.1 


14.9 


3.9 


3.1 


4.4 


1.688 


2 


17.2 


12.5 


20.2 


4.4 


3.4 


5.0 


1.194 


3 


20.5 


15.0 


24.1 


4.6 


3.5 


5.2 


0.844 


3 


22.0 


16.1 


25.7 


1.9 


3.8 


5.5 


< 0.8, n M = 


0.3 and £2a 


= 0.7 












6.754 





3.8 


3.3 


4.4 


2.0 


1.8 


2.2 


4.776 


1 


7.8 


6.1 


9.0 


2.7 


2.2 


3.0 


3.377 


2 


12.7 


9.6 


14.9 


3.4 


2.7 


3.9 


2.388 


3 


17.4 


12.9 


20.5 


1.0 


3.1 


1.6 


1.688 


3 


20.6 


15.2 


24.2 


1.6 


3.6 


5.2 


1.194 


3 


21.8 


16.1 


25.6 


1.8 


3.8 


5.5 



^Values of Lx min (0.5—2.0 keV) are based on the lower limits of the luminosity bins used in the derivation of the non-parametric XLF 



N mp = J X ' maX f ^ 0(Lx, z) n(f x {L x , z)) dz 

(7) 

Note that the non-evolving XLF is strictly a function of 
luminosity as parameterized by the Schechter fits to lo- 
cal clusters (e.g., Equation [S] with the best-fit parameters 
from Table However, we explicitly indicate the po- 
tential redshift dependence in this equation to generalize 
it for subsequent treatment with an evolving XLF. The 
statistical significance of the difference between iV b s and 
A^, xn is compu ted based on Poisson confidence intervals 
IGehreldllflsfo . 

Of the three local XLFs (BCS, RASS1BS, and RE- 
FLEX), we use the REFLEX measure as the preferred 
reference for three reasons: 1) the REFLEX local XLF is 
based on the largest sample of clusters used to date (452 
clusters), 2) the REFLEX normalization lies intermedi- 
ate to the BCS and RASS1BS, and 3) at low to inter- 
mediate luminosities, our 160SD low-redshift data most 
closely match the REFLEX normalization. Nonetheless, 
we will also quote our results relative to the BCS and 
RASS1BS to demonstrate the full range of possible sig- 
nificances. 

First we consider the redshift interval 0.3 < z < 0.6 
starting with the highest luminosity bin and then in- 
tegrating to lower luminosities in an Einstein-de-Sittcr 
model. In and above the highest bin of the XLF (Lx > 
4.8), we observe zero clusters whereas the REFLEX local 



XLF predicts 4.5 clusters according to Equation This 
difference is only 2.3a significant. The next two bins 
(L x > 3.4,2.4) have AT obs = 2,7 and iV cxp = 9.1,16.2 
with the associated significances of 2.6a and 2.4a. Con- 
tinuing from here to lower luminosity bins decreases the 
significance. These results are summarized in Tabled 
Examining the run of significance versus minimum lu- 
minosity without the constraints of the arbitrary bins 
of Figure El indicates that the significance briefly peaks 
above 3er near 3.4 x 10 44 ergs s _1 . However with the 
BCS XLF as the no-evolution reference, the apparent 
cluster deficit is not significant (i.e. always < 3d). 
Conversely, the RASS1BS XLF indicates a statistically 
strong signal across much of the probed luminosity. Re- 
peating this analysis in the context of a A-dominated 
universe uniformly reduces the significance of the clus- 
ter deficit such that only relative to the RASS1BS does 
the deviation seen in the 160SD measure at 0.3 < z < 0.6 
appear marginally significant (see Table Ej- 

Our highest redshift measure of the XLF, 
0.6 < z < 0.8, is clearly lower than all the determi- 
nations of the local XLF (Figures & [7J . Integrating 
over the entire luminosity range sampled in this redshift 
shell, a non-evolving model of the cluster population 
based on the REFLEX XLF predicts 22.1 clusters. This 
strongly conflicts with the actual observed sample of 
3 clusters, a 4.9a difference. The situation in the A- 
cosmology is essentially the same; 21.8 clusters expected, 
a 4.8a difference. If we instead use the RASS1BS XLF 
as our baseline, the predicted count is greater than 26 
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clusters with the resulting cluster deficit being 5.5er. 
More importantly, if we adopt the most conservative 
position (i.e., attempting to minimize the deficit) and 
use the BCS XLF, we expect to find 16 clusters which 
is 3.8a away from the observed value. An integrated 
bin-by-bin analysis is reported in Tabic 

4.2. Maximum Likelihood Analysis 

A more general approach to quantifying evolution in 
the cluster XLF is to perform a maximum likelihood fit 
of an evolving Schechter function to the observed cluster 
distribution in luminosity and redshift. This approach 
makes maximal use of the data, is free from arbitrary 
binning, and is sensitive to pote ntial negative evolution 
jRosati et al.ll2000tlHenrvll200^ . 

We follow the prescription of iMarshall et alJ l|1983f) 
but generalize the treatment to account for uncertain- 
ties in the observations. The luminosity-rcdshift plane 
is uniformly parsed into extremely small intervals of size 
dLxdz. In each element we compute the expected num- 
ber of clusters with luminosity Lx and redshift z: 



X(L x ,z)dL x dz = (j>(Lx,z) CL(fx(Lx, z)) —j-^- dLxdz. 

(8) 

The model XLF, c6(Lx, z), is an evolving Schechter func- 
tion of the form given in Equation [S] The important 
modification here is to free the density normalization, 
cf>* , and characteristic luminosity, L x , of the luminosity 
function to evolve with redshift such that, 



b*{z) 



and 



L*x(z) 



T* 



1 + z 

1 + Zq 



1 + z 



1 



zo 



(9) 



(10) 



where A and B parameterize the evolution. The local 
values of the normalization (</>g) and characteristic lumi- 
nosity (L x ) are taken from the local XLF determina- 
tion which samples a median redshift of Zq. This red- 
shift baseline is luminosity dependent as a result of the 
flux-limited nature of the survey technique. For instance 
in the BCS sample, at L x = 0.1 (1, 10) x 10 44 ergs s -1 , 
the median sample redshift is z = 0.02 (0.06, 0.21). 
Note that in the present model the faint-end slope of the 
luminosity function (a) is fixed at the local value since 
there is little evidence for any change in this parameter 
as a function of redshift. 

The sampling of the luminosity-rcdshift plane is suffi- 
ciently fine that in each element, dLxdz, the number of 
observed clusters is either zero or one and the expected 
number of clusters is much smaller than unity. In this 
sparse sampling limit we can define a likelihood function 
C based on joint Poisson probabilities, 



This is the combined probability of observing exactly 
one cluster at each point (£x,i> zt) populated by a 
160SD cluster and observing exactly zero clusters ev- 
erywhere else (Lxj,Zj). Again, occupied elements in 
the luminosity-redshift plane are indexed by i, whereas 
empty elements are indexed by j. Transforming to the 
standard expression, S = — 21n£, and dropping terms 
that are not model dependent, we have 



S = -2^w, ; ln[X(Lx,i,Zi)] 

\{L x ,z)dL x dz. (12) 



A weighting term, wt, is introduced to incorporate the 
uncertainties in observed cluster luminosity and the sum- 
mation is evaluated over a total number of i indices 
much greater than t he number of observed clusters (e.g., 
iBorgani et alJ 120011) . Instead of being a point in the 
observed (Lx, z)-plane, each cluster is smoothed in the 
Lx-direction according to a Gaussian distribution with 
a width set by the lcr flux error, ej, x , for each cluster 
(median value is 20% for the 160SD sample). No similar 
treatment is required for the redshift measures since the 
typical error is a few tenths of a percent. Thus a weight is 
assigned to each element in the luminosity-redshift plane 
based on the fractional contributions of all clusters in the 
same redshift interval, 



IV; 



E 



2-L 



exp 



{Lx,k - Lx^f 
2e?_ . 



dLy 



(13) 

Here the summation k is over the clusters with a redshift 
between zi — dz/2 and z\ + dz/2. 

Recall the overall goal is to find the values of A and B 
that predict a luminosity-redshift distribution that best 
matches the data. These best-fit parameters are deter- 
mined by minimizing S with confidence levels defined to 
be 

AS = S{A,B)-S(A hcst ,B hcst ) (14) 

Since S is distributed like x 2 , the lcr, 2a, and 3cr (68.3%, 
95.4%, and 99.7%) confidence intervals for a two p aram- 
eter fit are AS = 2.30, 6.17, and 11.8, respectively l)Avnil 
IT^IC^IT^IT97qh . 

We use the maximum-likelihood procedure to estimate 
the evolutionary parameters for a model XLF that best 
reproduces the observed cluster luminosity-redshift dis- 
tribution. For the local parameters of the XLF Lx,o> 
and a) we adopt the values from the BCS. This per- 
mits comparisons to previous work and assumes a con- 
servative position since the BCS has the lowest normal- 
ization of the three local references. We estimated the 
baseline redshift, zq(Lx), by calculating the median red- 
shifts in seven luminosity intervals of the BCS sample 
and interpolating across these points (zq = 0.018, 0.028, 
0.040, 0.061, 0.103, 0.171, 0.236 at L x ,Eds/10 44 = 0.054, 
0.139, 0.303, 0.938, 2.283, 5.771, 12.929 ergs s" 1 ). The 
top panels of Figure [S] indicate the constraints using the 
entire 160SD sample of 177 clusters (0.02 < z < 0.80, 
10 42 < L x < 10 46 , fx > 4 x 10~ 14 ) in two dif- 
ferent cosmologies: A = 2.6±g;?, B = -3.2±°:e (EdS), 
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Fig. 8. — Maximum-likelihood contours for the parameters A and B characterizing the evolution of the cluster XLF where </>* <x 0q(1 + z) A 
and <x Z/J (1 + z) B based on four independent surveys. Plotted are the results from our analysis of the 160SD and NEP samples, plus 

similar data for the EMSS and RDCS from lRosati et alj |2002). The no-evolution scenario corresponds to A = B = 0. Confidence regions 
in the left column are computed in an Einstein— de-Sitter universe, whereas those in the right column are for a A-dominated universe. TOP 
PANELS: One, two, and three a contours with the 160SD results based on the entire sample. MIDDLE PANELS: One, two, and three <r 
contours with the 160SD results based on the z > 0.3 sample. BOTTOM PANELS: Same as the middle panels except only la contours 
are shown for clarity. 



A = 1.3±° 6 7 , B = -2.3±° l (A-cosmo). In either case, no 
evolution (A = B = 0) is strongly excluded at > 3a. 

To examine the evolution "signal" in the 160SD sam- 
ple at higher redshifts corresponding to the XLFs in 
Figures & we repeat the maximum-likelihood anal- 
ysis using the 66 clusters at z > 0.3. The resulting con- 
tours are shown in the middle panels of Figure |5] and 
the best-fit values are: A = 1.6± ;g, B = -2.9±^ (EdS), 

A = 0.6+?;5\ B = -2A± X q% (A-cosmo). Relative to the 



full sample analysis, the best-fit A is about la lower 
while the best-fit B is effectively unchanged, being a few 
tenths of a lower. Apparently, the A parameter in the 
full analysis is elevated to fit an enhancement of clusters 
at z < 0.3 which can be seen at ~ 3x 10 43 ergs s _1 in 
the low-redshift XLF. We adopt the less-biased z > 0.3 
results as representative of the 160SD in subsequent dis- 
cussions. 

For comparison purposes, we also show in Figure [Slthc 
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Fig. 9. — Comoving volume density of very X-ray luminous 
clusters with Lx > 10 45 ergs s _1 (0.1— 2.4keV; equivalently 
>6.2 X 10 44 , 0.5-2.0keV). The data points are derived from the 
eBCS and t he MACS samples as reported by H. Ebeling et al. (see 
Figure 5 in Ellis 2002, and H. Ebeling et al. 2004, in preparation) 
in an Einstein-de-Sitter universe. The dashed line and shaded re- 
gion indicate the prediction based on our ensemble best-fit model 
XLF (A = 0, B = -2.5 ± 1); not a fit to the eBCS+MACS data. 
The dotted line marks the no-evolution expectation derived from 
the BCS Schechter function fit to the local XLF which is dominated 
by z ~ 0.2 clusters at these luminosities. 



constraints from three additional datasets — the NEP, 
EMSS, and RDCS. One, two and three a contours are 
shown in the top and middle panels, whereas only la 
contours are plotted in bottom panels. The degree of 
agreement amongst the four surveys is reasonably good 
but certainly not perfect. It is not immediately clear 
how much of this is due to potential systematics in the 
individual datasets or the appropriateness of the model 
XLF used in the maximum-likelihood fit. These issues 
will be addressed in a forthcoming paper focused on the 
joint analysis of the available samples (J. P. Henry et 
al. 2004, in preparation; preliminary results discussed in 
IHenrvll2003T) . 

Considering the ensemble of data from four in- 
dependent cluster samples, the no-evolution scenario 
(A = B = 0) is strongly ruled out at 3> 3a at z < 0.8 re- 
gardless of cosmology. In an Einstein-de-Sitter universe, 
the change in the cluster population is consistent with 
pure luminosity evolution (A rj 0,5 rj —2.5), whereas 
it may be a combination of luminosity and density evo- 
lution (A w — 1,5 « —2) in a A-dominated universe. 
Given the lack of clusters with Lx > 10 45 ergs s _1 in 
the analyzed samples, caution is warranted if these re- 
sults are extrapolated to predict the cluster abundance at 
very high luminosities. Nonetheless, we show in Figure 
that our characterization of cluster evolution is consis- 
tent with a direct measurement of the comoving volume 
density of clusters at L x > 10 45 from the eBCS+MACS 
surveys (H. Ebeling et al. 2004, in preparation). 

5. CONCLUSIONS 

We have used the 160SD sample of 201 X-ray selected 
galaxy clusters to track the volume density of these sys- 
tems at local, intermediate and high rcdshifts. Nonpara- 
metric measurements of the cluster XLF suggest there 
is effectively no detectable evolution in the population 
out to z ~ 0.6 at luminosities less than a few times 



10 44 h^Q ergs s _1 . However, data in the redshift inter- 
val 0.6 < z < 0.8 indicate a mild but significant (> 3a) 
change in volume densities above luminosities of approx- 
imately 10 44 h^Q ergs s _1 . Our findings demonstrate a 
deficit of high-luminosity clusters at high redshift rel- 
ative to the present-day levels. For example, we ob- 
serve only 3 clusters at 0.6 < z < 0.8 where an integral 
over the local XLF predicts we should find at the very 
least 16 clusters. A maximum-likelihood analysis of the 
observed luminosity-redshift distribution further under- 
scores that the no-evolution scenario is entirely inconsis- 
tent with 160SD data. Modelling the XLF with an evolv- 
ing Schechter function, we have demonstrated that the 
160SD, NEP, RDCS and EMSS clusters samples all re- 
ject a model lacking evolution. Further our evolving XLF 
model derived from datasets at Lx < 10 45 successfully 
predicts the observed volume densities at Lx > 10 45 . 

A composite view of the high-redshift cluster popu- 
lation is shown in Figure ITUI Here we plot the latest 
compilation of distant XLFs as measured by eight X-ray 
flux-limited surveys. From these data six independent 
teams have co ncluded there is mea s urable cluste r evolu - 
tion — EM SS:lc!oTaet~aT| lll990alL IHenrv et all lf l992L 
BSHA RC: INichol et al.l~ i 1999T . RDCS: IRosati etall 
( 2000), NEP: lGioia et all (|2001f) . 160SD: this work, and 
MACS (H . Ebeling 2004, in preparation). Conversely, 
SSHA RC HBurke et al.M1997t) and WARPS Hones et alJ 
2000a) draw the opposite conclusion. Note however that 
the sky coverage of the SSHARC survey (17.7 deg 2 ) is 
inadequate to measure evolutionary changes at the very 
bright end of the XLF, as studied in the large solid angle 
EMSS (735 deg 2 ) and the 160SD surveys. 

Measurements of the cluster X-ray luminosity function 
are consistent with a population whose comoving vol- 
ume density evolves as a function of both redshift and 
luminosity (mass). For a sample of distant clusters in a 
fixed luminosity interval, the difference in volume den- 
sity for this sample relative to the present-day popula- 
tion increases with increasing redshift (e.g., Figure 0). 
Considering a sample of clusters in a distant redshift in- 
terval, the difference in volume density for this sample 
and the present-day population increases with increas- 
ing luminosity (e.g., Figure ITU)) . Visualizing this on the 
luminosity-redshift plane, the cluster deficit is maximal 
at high redshift and high luminosity (i.e., the upper-right 
corner of Figure 0. The maximum-likelihood results 
in Figure El demonstrate that the EMSS, 160SD, RDCS, 
and NEP surveys all reject the no-evolution scenario and 
generally agree on the character of the evolution (e.g., 
AB contours overlap at 1-2 a). The EMSS was able 
to detect a statistically significant cluster deficit in the 
0.3 < z < 0.6 redshift shell because of their sensitivity to 
very luminous clusters. Surveys such as the 160SD probe 
lower luminosities where the evolution "signal" is intrin- 
sically weaker and thus have to push to higher redshifts 
(e.g., 0.6 < z < 0.8) before the deficit grows large enough 
for a significant detection. 

The preponderance of observational data requires mild 
evolution in the bright end of the cluster XLF as mea- 
sured out to z ~ 0.8, whereas the volume densities of the 
bulk of the population (low and intermediately massive 
systems) do not change appreciably. In the context of hi- 
erarchical structure formation, we are probing sufficiently 
large mass aggregations at sufficiently early times where 
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Fig. 10. — Compilation of high-redshift XLFs as measured by eight i ndependent X-ray fl ux-limi ted surveys fEinst ein-de-Sitt er universe). 
Refere nces are 160 SD (this paper), R DCS IRosati et aLll^OOOT). EMS S IHenrv et ai]IT99l) . NEP IGioia et aljl200ll) . WARPS Uones et all 
BOQOhfl . SSHARC IBurke et nUMm . BSHARC INichol et alllT99Sh . and MACS (H. Ebeling 2003, private communication). The shaded 
region delineates the regime of the local XLF whereas the heavy solid line is an evolving model XLF (A = 0, B = —2.5) at z = 0.7. 



the Universe has yet to assemble these clusters to present- 
day volume densities. CosmoloKical constraints derived 
from the RDCS and 1 60SD samples l|Borgani et all200H 
Vikhl mTn et alJ l2003|) demonstrate the observed evolu- 
tion is consistent with structure growth in a flat universe 
with J2m ~ 0.3. 

A decade after the first evidence for negative clus- 
ter evolution reported by the Einst ein EMSS survey 
l|Gioia et alJll990at IHenrv et all ll992). we are now wit- 
nessing the fruition of the ROSAT era cluster surveys: 1) 
the robust determination of the local cluster X-ray lumi- 
nosity function, and 2) multiple, independent confirma- 
tions of the EMSS results on cluster evolution. Further 
work remains to assimilate the available cluster samples 
in a detailed, joint analysis to produce maximal con- 
straints on the evolution phenomenon. Furthermore, 
additional inputs are expected in the near term from 
two on-going surveys — BMW-HR I ijPanzera et all2003f) 
and MACS labeling et al.ll20f)l . The BMW-HRI has 
a sky coverage rather similar to the 160SD, and hence 
will likely provide similar XLF measures. The final re- 



sults from MACS should prove particularly interesting 
because this program probes extremely high luminosi- 
ties (Lx ^ 10 45 ergs s _1 ) at high redshift exactly where 
the evolution signature should be strongest. Finally, on 
longer timescales, searches with XMM-Newton and pos- 
sibly a new dedicated X-ray survey satellite should reveal 
substantial numbers of z > 1 clusters which will provide 
powerful leverage in the study of cluster evolution. 
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